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ABSTRACT 

This paper summarizes a research program that has been underway for a decade. The 
objective is to find a fast and accurate scheme for solving quantum problems which does 
not involve a Monte Carlo algorithm. We use an alternative strategy based on the method 
of finite elements. We are able to formulate fully consistent quantum-mechanical systems 
directly on a lattice in terms of operator difference equations. One advantage of this 
discretized formulation of quantum mechanics is that the ambiguities associated with op- 
erator ordering are eliminated. Furthermore, the scheme provides an easy way in which to 
obtain the energy levels of the theory numerically. A generalized version of this discretiza- 
tion scheme can be applied to quantum field theory problems. The difficulties normally 
associated with fermion doubling are eliminated. Also, one can incorporate local gauge 
invariance in the finite-element formulation. Results for some field theory models are 
summarized. In particular, we review the calculation of the anomaly in two-dimensional 
quantum electrodynamics (the Schwinger model). Finally, we discuss nonabelian gauge 
theories. 

PACS numbers: 02.90.+p, 11.90.+t, ll.lO.-z 
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I. INTRODUCTION 



A classical physical system is described by a differential equation supplemented by ini- 
tial conditions. The differential equation characterizes the time evolution of the system. In 
quantum physics one cannot specify initial conditions for an operator differential equation 
because specifying the values of the fundamental observables (operators) at a given time 
would violate the uncertainty principle. Thus, a quantum theory is described by an oper- 
ator differential equation (a field equation) supplemented by an equal-time commutation 
relation (ETCR). The ETCR is a constraint which must hold at each time; its constancy 
in time expresses the unitarity (conservation of probability) of the quantum theory. The 
ETCR of quantum mechanics replaces the initial condition of classical mechanics. 

In the study of quantum field theory there are inherent ambiguities and divergences 
associated with an operator-differential-equation formulation because operator-valued dis- 
tributions in the continuum are so singular that products of such operators do not exist. It 
is well known that introducing a space-time lattice is a good way to remove these ambigu- 
ities and thus to regularize a continuum quantum field theory. The content of the theory 
is then contained in the continuum limit of the lattice theory. 

Ordinarily, the lattice is introduced as a mathematical artifice to make sense of the 
functional integral representation of a quantum field-theory. In this computational scheme 
the lattice regularizes the functional integral as an infinite product of ordinary Riemann 
integrals. Then, the infinite product is approximated as a finite product of integrals which 
are evaluated by Monte-Carlo methods. This procedure is slow; doubling the computer 
time gives only minimal improvements in accuracy. 

The program discussed here uses the lattice in a completely different and more funda- 
mental way. We show how to formulate and construct a fully consistent (unitary) quantum 
theory on a space-time lattice. Such a theory is defined in terms of an operator difference 
equation (rather than a differential equation) and an ETCR that holds at equal-time lat- 
tice points. The operators in such a theory have none of the problems (infinities) that 
operators in the continuum have. We will see that, while there is no hope of solving op- 
erator differential equations, operator difference equations can be solved exactly. Thus, 
for each lattice we obtain exact closed-form solutions for the field operators rather than a 
slowly converging sequence of Monte-Carlo approximations. 

To convert an operator differential equation to an operator difference equation we use 
the method of linear finite elements. This method is explained in Sec. II of this paper. 
There we show explicitly that for a one-particle quantum system the discrete-time operator 
equation obeys the constraint of unitarity. We show how to solve the operator difference 
equation exactly and how to use the solution to obtain accurate numerical estimates of 
the eigenvalues. In Sec. Ill we show that the method of finite elements resolves the well- 
known operator-ordering ambiguities one encounters in the usual Hamiltonian formulation. 
In Sec. IV we show that one can use higher-order finite elements to generate systems of 
operator difference equations. We show that the requirement of unitarity can be used 
to recover gaussian quadrature. In Sec. V we apply the method of finite elements to 
Hamiltonian systems, such as spin systems, that are associated with algebras other than 
the Heisenberg algebra. 
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The next five sections address the more difficult problem of systems having more than 
one degree of freedom. In Sec. VI we consider the case of two degrees of freedom and 
show that on a finite-element lattice the discrete quantum system is unitary. Next, we 
examine the question of how to apply the method of finite elements to quantum field 
theory problems, which are systems having an infinite number of degrees of freedom. In 
Sec. VII we consider the simplest case; namely, that of scalar quantum field theory. Then, 
in Sec. VIII we examine the Dirac equation and spinor field theories. We show that 
the method of finite elements resolves the fermion-doubling problem, a generic difficulty 
usually encountered whenever one attempts to discretize the Dirac equation. In the last 
two sections we use the method of finite elements to solve quantum field theories exhibiting 
local gauge invariance. In Sec. IX we examine quantum electrodynamics, a quantum field 
theory that has an abelian gauge invariance. We apply our analysis to solve the Schwinger 
model, massless quantum electrodynamics in two dimensions, and obtain the anomaly. In 
Sec. X we examine quantum field theories having local nonabelian gauge invariance. 
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II. QUANTUM MECHANICS WITH ONE DEGREE OF FREEDOM 



A. The Method of Finite Elements 

The method of finite elements^ is a technique for solving partial differential equations 
that is well known to applied mathematicians. The method consists of four steps. We are 
given a classical partial differential equation Lcj) = 0, which is to be solved on a region R 
subject to boundary conditions given on the boundary OR. We first decompose R into 
a set of nonoverlapping patches, called finite elements, which completely cover R. For 
classical (not quantum) differential equations the patches may have arbitrary sizes and 
shapes. On each patch we approximate the solution to the partial differential equation 
as a polynomial. The degree of this polynomial is chosen to suit the conditions of the 
problem. Second, at the boundaries of contiguous patches continuity is imposed (and 
sometimes continuity of higher derivatives). Third, on patches that are adjacent to dR we 
impose the boundary conditions. Fourth, we impose the differential equation Lcp = at 
one point (or more than one) on each patch. Conditions two, three, and four give a system 
of algebraic equations satisfied by the coefficients of the polynomials. Solving this system 
gives a good approximation to the solution 0. 

We illustrate this procedure by solving a simple classical ordinary differential equation 
problem: 

y\x)^y{x), 2/(0) = 1. (2.1) 

Show that 

y(l) =e = 2.71828.... (2.2) 

We begin by using just one linear finite element: y = ax + b, where < a; < 1. The 
initial condition gives one algebraic equation. 

y(0) = l 6=1. (2.3) 

We must impose the differential equation at one point xq on the interval < a;o < 1; 
however, the choice of xq remains ambiguous. Later we will see that unitarity in quantum 
mechanics removes this ambiguity and uniquely selects xq = 1/2. For now we simply 
choose a;o = 1/2 and proceed: 

2/'(l/2) = 2/(1/2) ^ «=i + ^- (2-4) 

Solving (2.3) and (2.4) for a and b gives y{x) = 2a; + 1, so that y{l) = 3, which is a good 
result that already differs from the exact answer by only 10%. 

For the case of two linear finite elements yi — at + b, y2 — ct + d, where t is a local 
variable that ranges from to 1/2, the initial condition gives 

2/i(0) = l ^ b=l. (2.5) 

Continuity at x = 1/2 gives 

yi(l/2) = 2/2(0) ^ + b = d (2.6) 
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and imposing the differential equation at the center of each finite element x = 1/4, and 
a; = 3/4 gives 

yi(l/4) = yi(l/4) ^ «=^ + ^ 

2/2(1/4) = y2 (1/4) ^ c=| + d. (2.7) 
Simultaneous solution of Eqns. (2.5-7) gives an excellent result for y(l): 

2/(1) =2/2(1/2) = y = 2.778..., (2.8) 

which differs from the correct answer by 2%. 

In general, for N finite elements the exact result for the approximate value of y{l) is 

.(i)=(^r, (2.9) 



2N -1^ 

which for large N approximates e with a very small relative error of 1/(12A^^). 

B. The Equations of Quantum Mechanics 

Applying the technique of finite elements to quantum problems is much more inter- 
esting because the polynomial coefficients are operators. Consider the simple quantum- 
mechanical Hamiltonian 

H=\p' + y{(i). (2.10) 

for which Hamilton's equations are 

q = P. P=-V'{q). (2.11) 

The system (2.11) constitutes a time-evolution problem for the operators p{t) and q{t). 
The analog of the classical initial condition is an operator constraint in the form of an 
ETCR 

[q{t),p{t)] = i. (2.12) 

If (2.12) is imposed at t = then, by virtue of (2.11), it holds for all t. 

To solve (2.11) on the interval [0,T] we introduce a lattice of N linear finite elements. 
On each finite element t ranges from to h and Nh = T. Let g^, ("- = 0, 1, 2, . . . , A^) 
be the approximate value of q{nh). Let us examine the nth finite element, where p{t) and 
q{t) are approximated by the linear polynomials 

p{t) = (1 - t/h)pn-l + it/h)pn, 



q{t) = (1 - t/h)qn-i + {t/h)qn, (2.13) 
5 



with < t < h. Substituting (2.13) into (2.11) and evaluating at the center of the finite 
element t = h/2, we obtain a pair of algebraic equations relating the operators Pn-i, Qn-ij 
Pn, and Qn. 

Qn - Qn-l Pn + Pn-1 in ^ A \ 

-nr- = -^—^ ^^-'^^^ 

Pn - Pn-1 ^ _y, f gn+gn-A ^2.14b) 



h V 2 ^ 

The ETCR (2.12) at the lattice point nh reads 

[Qn.Pn] = i- (2.15) 

It is not obvious that (2.14) and (2.15) are consistent. To prove consistency we argue as 
follows: Equation (2.14a) implies that 

[qn-qn-l,Pn+Pn-l]^0 (2.16a) 

and (2.14b) implies that 

[qn+qn-l,Pn-Pn-l]=Q- (2.16b) 

Expanding and adding together the commutators in (2.16) gives the result 

[Qu.Pu] - [qn-UPn-\\ = 0. (2.17) 

Thus, if [qojPo] = i initially, then (2.15) holds for all values of n. The proof of the 
persistence of the ETCR's holds if and only if both differential equations (2.11) are imposed 
on finite elements t — h/2. At every other point on the finite element (2.15) ceases to be 
true.^ Thus, quantum-mechanical unitarity (persistence of the ETCR's) removes a basic 
ambiguity that occurs in the numerical solution of classical differential equations; namely, 
where on the finite element to impose the differential equation. 

C. Solution of the Operator Equations 

We have proved the consistency of (2.14), but we must now solve these difference 
equations. To do so we use (2.14a) to eliminate pn from (2.14b) which now becomes 

4 ^2 ^.,fqn + qn-i\ , 4 fqr, + qn-i\ .^,„s 

If we let X = {qn + qn-i)/2, y = Aqn-i/h? + 2pn-i/h, and g{x) = V'{x) + Ax/h"^, then 
(2.18) becomes 

y = g{x). (2.19) 

While X and y are operators, (2.19) implies that they commute and thus (2.19) can be 
treated as if it were a c-number equation. Its exact solution is a; = g~^{y) so 



-1/^2 4 \ 

qn = -qn-l + 2g I J^Pn-l + j 



Pn = -Pn-1 



(2.20) 



This result shows that the exact operator solution after N time steps to the lattice quantum 
theory in (2.14-15) is a continued (nested) function.^ 

The unitarity of the lattice theory can also be demonstrated explicitly because the 
transfer (lattice time evolution) operator U can be expressed in closed form: 



qn+i = Uq„U \ Pn+i = Up„U \ 



where 



with 



U 



eihpl/4^ihA(qr,)^ihpl/4 



A{x) 



2 



+ V 



(2.21a) 



(2.21b) 



It is interesting that while the solution in (2.20) and the transfer operator U involves 
the function g~^, matrix elements of these operators only involve g. For example, if we 
define Fock states \n) at the initial time by 



Pq 



1 



(a — a^) and go 



^(a + at), 



where 



then/ 



an 



) — \/n\n — 1) and a^\n) = \/n-\- l|n+ 1) 



{m\qi\n) = ^(x/n5^,n_i + \/rn5n,m-i) 



V2 



gi6(n— m) 



i?V7r2'^+'^n! 
where B? = 47^/i~^ + /i~^7" 

D. Energy Eigenvalues 



2R 



9M 
2R 



cos 6* = 27/(i?/i^), and Hn is the nth Hermite polynomial. 



It is easy to compute energy levels of quantum systems once the operator equations 
have been solved. There are two ways to carry out such a calculation. The quick and 
approximate method makes use of the one-finite-element solution to the operator difference 
equations. The techniques for performing this calculation are given elsewhere.^ Here are 
some numerical results: For the case of the harmonic oscillator, the energy gap u = Ei — Eq 
comes out exactly. For the anharmonic oscillator, where V{q) = the exact value 

of uj is 

1.08845... AV3_ The one-linear-finite-element equations predict 1.14471A^/^ (5.2% 
relative error); the one-quadratic- finite-element equations (see Sec. IV) predict 1.08225A^/^ 
(-.57% relative error). 

It is also possible to determine all energy differences simultaneously and to arbitrary 
accuracy by taking large numbers of finite elements. The procedure consists of using the 
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method of finite elements to calculate the matrix element An = {0\qn\l) of the operator 
Qn- The sequence of numbers is a discrete time sequence. We can then compute the 
discrete Fourier transform A^, of this sequence and search for peaks in this distribution. 
These peaks correspond to energy differences in the spectrum of the theory. The detailed 
procedure for this calculation is given elsewhere.^ 

Figure 1 shows the results of a short computer calculation using 1000 finite elements 
for the anharmonic oscillator. The results in Fig. 1 are summarized numerically in Table 

1.6 



Figure 1. A semilog plot of versus m for the anharmonic oscillator V{q) = 0.885g^. 

The spikes give extremely accurate approximations to energy differences Ej — of the 
exact spectrum. To read off the predicted energy differences, we note that one unit on 
the horizontal scale corresponds to an energy increment of A = 2tt/[{N -\- l)h]. Energy 
differences are measured from both the left axis and the right boundary (see Ref. 6). 
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Energy difference Exact Approximate Reiative error 



El — Eq 


1.728 


1.674 


3.1% 


E2 — El 


2.142 


2.218 


-3.5% 


Es — E2 


2.537 


2.595 


-2.3% 


Ea — E3 


2.790 


2.846 


-2.0% 


E5 — E4 


3.000 


3.097 


-3.2% 


Eq — E^ 


3.210 


3.306 


-3.0% 


hjo — rL/Ci 


6.407 


6.487 


-1.2% 


E/i — El 


7.469 


7.659 


-2.5% 


E5 — E2 


8.327 


8.454 


-1.5% 


Eq — E3 


9.000 


9.165 


-1.8% 


E5 — Eq 


12.20 


12.39 


-1.6% 


Eq — El 


13.68 


13.98 


-2.2% 



Table 1. Comparison between exact eigenvalue differences for the anharmonic oscillator 
V{q) = gq^, with g — 0.885 and the approximate eigenvalue differences obtained in the 
following manner. We find the integer values of m for which |74(m)p is a local maximum. 
The energy difference Ej — Ef., for some j and k, is then predicted to be 27rm/[{N + 
l)h]. This procedure gives relative errors of order 1-3% as the table shows. However, 
this procedure can be drastically refined by using by using nearby values of |^(m)p to 
interpolate the precise noninteger value of the location of the maximum. 

Computations involving many finite elements also give extremely accurate results in 
problems involving tunneling.^ Rather than solving the time-dependent Schrodinger equa- 
tion we can obtain quantum tunneling results directly by following the time evolution of 
the matrix elements of operators on a discrete-time lattice. 
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III. GENERAL HAMILTONIAN SYSTEMS AND OPERATOR ORDERING 



The Hamiltonian H in (2.10) is not the most general one-particle quantum-mechanical 
Hamiltonian. However, the more general one-particle Hamiltonian, H{p, q), suffers from an 
ambiguity that is not present in H in (2.10); namely, the well-known problem of operator 
ordering. There is no such ordering problem in classical mechanics. However, quantizing a 
given classical-mechanical Hamiltonian H{p, q) is an ambiguous procedure. The following 
simple example illustrates why this is so. Consider the classical Hamiltonian H{p, q) = 
p^q^. The corresponding quantum- mechanical Hamiltonian could consist of any of the 
following Hermitian operators: 

^(Z^^ + ^V), ^iPQPQ + QPQP), PQ^P, QP'^Q- 

It could even be a linear combination of these operators. We show in this section that if we 
take the lattice as fundamental and use the method of finite elements, the above operator 
ordering ambiguity is completely eliminated.^ 

With H{p, q) an arbitrary function of the operators p and q the operator differential 
equations of motion in the continuum are 

Pit) = = -i\p, H]. (3.1) 

Let us address the problem of converting this system of equations in the continuum 
into a system of unitary operator difference equations on a time lattice. Recall that by 
unitary we mean that the difference equations exactly preserve the equal-time commutation 
relation 

m.pit)]=i (3.2) 

at each time step. We show that if the method of finite elements is used to construct 
the operator difference equations then the ordering of the operators p and q is uniquely 
determined by the unitarity requirement. 

The finite-element method consists of making the replacements 

qit) (qn - qn-i)/h, 

pit) {pn-Pn-l)/h, 
qit) (<?n +gn-l)/2, 

p{t)^{Pn+Pn-l)/'2, (3.3) 

in (3.1). Thus, on the lattice, the differential equations (3.1) become 

Q^^H{P,Q), P = -±H{P,Q), (3.4) 
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where we have used the notation 



Q = {qn- qn-i)/h, 

P = {Pn -Pn-l)/h, 
Q={qn + Qn-l)/2, 

P=(p,+p,_i)/2. (3.5) 
To estabhsh unitarity as in Sec. II one must prove that 

[Qn,Pn\ = [qn-l,Pn-l], ("-=0, 1, 2, . . .). (3.6) 

We can establish (3.6) if we can explicitly show from the operator difference equations 
(3.4) that 

[Q,P] + [Q,P]=0. (3.7) 

To see this we substitute the definitions in (3.5) into (3.7) and expand the commutators; 
this calculation directly shows that (3.7) implies (3.6). Thus, our objective is to examine 
the expression 

[-^H{P, Q),P]- [Q, ^H{P, Q)] (3.8) 

and to show that it vanishes. 

It is crucial to remark that in general (3.8) does not vanish. This is because the 
commutator 

e=[Q,P] (3.9) 

is not a c- number even though [q(t),p{t)] is; the quantity 9 is an operator because it 
contains the unequai-time commutators [qn-i,Pn] and [qn,Pn-i]- 

To investigate (3.8) we may assume that H is Hermitian and that H{P, Q) can be 
expanded in a series of Hermitian terms Hm,n{P, Q): ip^, n >0), which consists of a sum 
of monomials containing m factors of P and n factors of Q. We can examine each term 
Hm,n of the series independently. For example, H2,2 has the form 

H2,2 = aPQ^P + bQP^Q + c{P^Q^ + Q^P^) + d{PQPQ + QPQP), (3.10) 

where a, b, c, and d are real constants. To illustrate our procedure we examine i?2,2 in 
detail. We compute 

= (2c - a - d){dQP + PQ0) + (a - b){Q9P + PdQ) + {b + d- 2c){0PQ + QP0). (3.11) 

Thus, unitarity requires that 2c — a — d= 0, a — b = 0, and b + d — 2c= 0. The solution to 
these equations is a = 6 and d = 2c — b where b and c are arbitrary real constants. Thus, 
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it appears that there is a two-parameter family of Hamiltonians of the type H2^2 which 
exhibit unitarity on the lattice. Indeed, H2.2 can be written in the form 



H2,2 = CT2,2 + (C - 6)G2,2, (3.12) 

where 

T2,2{P, Q) = PQ^P + QP^Q + P^Q'' + Q^P"" + PQPQ + QPQP (3.13) 

and 

^2,2 = PQPQ + QPQP - PQ^P - QP^Q. (3.14) 

Thus, if we were given a continuum Hamiltonian of the general form if2,2(p, q) we 
could reorder the operators p and q [using the commutation relation (3.2)] to make it take 
the form in (3.12) before making the finite-element transcription (3.3). [Of course, this 
reordering of operators produces additional simpler terms of the Hi i{p, q) = a{pq + qp)-] 
However, we observe that G2,2{p, q) is trivial; using (3.2) we see that ^2,2(^5 q) = —1- The 
above calculation shows that the requirement of lattice unitarity forces us to preorganize 
the operators p and q in H2^2{p,q) in a unique way; namely, the totally symmetric sum 
(T form) in (3.13). For example, if we are given the Hamiltonian i/2,2 = bqp^q, this 
Hamiltonian must be (uniquely) reordered by using (3.2) as 

H2,2{p,Q) = It2,2{p,<i) + 1 (3-15) 

before going onto the lattice by use of (3.3). 

What is remarkable is that given any Hamiltonian H^^^iP, q) carrying out the above 
procedure shows that there is always a unique form which is necessary and sufficient in 
order that the equal-time commutators be preserved as in (3.6). In particular, we must 
rewrite 

Hm,n{P^ 1) = 0:Tm,n{P, Q) + Hm-2,n-2{p, q), (3.16) 

where Tm,n is the totally symmetric sum (T form) of all possible monomials containing m 
factors of p and n factors of q. This process is then iterated until Hm,n is a descending 
sum of totally symmetric parts: 

Hjn,n{p,q) = aTm,n{p,q) + PTjn-2,n-2{p,q) + ■ ■ ■ ■ (3.17) 

To verify this assertion we must use the fact that derivatives leave the T form intact. 
In fact we have the identities 

-^Tm,n{P^ Q)={m + n)Tm-l,n{P^ Q), 

^Tm,n{P,Q) = {m + n)Tm,n-i{P,Q). (3.18) 
In addition, we observe that commutators maintain the totally symmetric form 

[Q,Tm,n{P,Q)] = Tm-l,n,l{P,Q,0), 
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[T^,„(P,Q),P] = T^,,_i,i(P,Q,^), (3.19) 

where Tm,n,i{P-i Q-i ^) is the totally symmetric sum of all monomials having m factors of 
P, n factors of Q, and one factor of 9. Using (3.18) and (3.19) it is easy to verify that the 
expression in (3.8) vanishes when H{P, Q) is in T form. 

This ordering procedure applies to all Hamiltonians H{p, q) which are polynomials in 
the variables p and q. However, if is a nonpolynomial function the ordering problem is 
much more challenging. For example, consider a class of Hamiltonians of the form 

H{p, q) = H{T,,,) = H{pq + qp). (3.20) 

To order the operators of this Hamiltonian we introduce a little-known set of orthonormal 
polynomials Sn{x) called continuous Hahn polynomials.^'^°'^^ These polynomials emerge 
from the simple observation that T„^^ is a polynomial function of Ti i; the defining equation 
for Sn{x) is therefore^^'^^ 

-5n(Ti,i) = r„,„/(2n-l)!!. (3.21) 
The first few polynomials Sn{x) are 

Sq{x) = 1, 
Si{x) = X, 

S2{X) = 1{X'-1), 

1 „ 

Ssix) = -{x - 5x), 
S^{x) = ^{x^-Ux^ + 9), 
S^ix) = ^(x^-30x3 + 89x), 

Sq{x) = ^{x^ - 55x^ + 439a;^ - 225). 

These polynomials have the following properties: 

(i) The generating function G{t) is 

„arctant 

^^^^^ (l + mi/2 -T.Sr.{x)t\ (3.22) 

^ ' n=0 

(ii) The orthonormality condition is 

dx w{x)Sm{x)Sn{x) = 5mn, (3.23) 
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where the weight function w{x) is given by^^ 

«;(a;) = [2cosh(7ra;/2)]-^ (3.24) 

(iii) A recursion relation satisfied by Sn{x) is 

nSn{x) = xSn-i{x) - (n - l)Sn-2{x). (3.25) 

The polynomials do not satisfy a second-order differential equation but they do obey the 
second-order functional difi'erence eigenvalue equation 

(1 - ix)Sn{x + 2i) + (1 + ix)Sn{x - 2i) = (4n + 2)Sn{x). (3.26) 

Now we return to the problem of ordering the operator H in (3.20). Using the com- 
pleteness of Sn{x) we expand H{x) as a series in Sn{x): 

oo 

H{x) = Y,(^nSn{x), (3.27) 



n=0 



where 



Thus, from (3.21) we have 



Q'n = I dx w{x)H{x)Sn{x). (3.28) 



H{Ti,i) = Yl ^nTn,n/{2n - 1)!!. (3.29) 



n=0 



We have therefore represented H{pq + qp) as an infinite sum of operators in T form. In 
the form (3.29) H is directly amenable to lattice transcription and the resulting operator 
difference equations automatically preserve unitarity. 

As an example, consider the Hamiltonian H = e'^^^'^~^'^^\ where c is a constant. For the 
exponential function, the integral in (3.28) can be performed in closed form and the result 
is an = (tanc)''[l + (tanc)2]i/2. Thus, 

oo 

H=[l + (tanc)2]i/2 ^(tanc)"5n(pg + qp) 

n=0 

= [l + (tanc)^]V2f;i^5^T„^^. (3.30) 

n=0 ^ ^" 

H is now in its unique T form and therefore the resulting Heisenberg equations can be 
transcribed onto the lattice. 
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IV. HIGHER-ORDER FINITE ELEMENTS 



It is possible to generalize the linear polynomials representing p{t) and q{t) to polyno- 
mials in arbitrary degree r: 



p{t) = J2Mt/h)\ (4.1a) 



k=0 

r 

q{t) = Y,hk{t/h)\ (4.1b) 

It is now necessary to determine r + 1 pairs of coefficients on each finite element interval. 
The procedure for doing so is evidently ambiguous; if we impose the differential equations 
(2.11) d times on the nth interval, then it is necessary to impose r-|- 1 —d joining conditions 
(continuity, continuity of the first derivative, continuity of the second derivative, and so 
on) at t = {n — l)h. On the first interval, there are no joining conditions at t = 0; 
rather we must impose r -\- 1 — d initial conditions, in which the values of q'(O) and 7j(0), 
q'(O) and p(0), g(0) p(0), and so on, are specified. These values are obtained by 
successively differentiating the differential equations (2.11). We say that as the number of 
joining conditions increases the approximation becomes stiffer. In one extreme, the stiffest 
approximation, the differential equation is imposed once on the interval, and in the other 
extreme, the floppy approximation, the method we will use in this paper, the differential 
equation is imposed r times, and we require only that the approximation be continuous. 

A. Failure of the Stiff Approximation 

The stiff approximation is forbidden by quantum mechanics. For a quantum-mechanical 
system with operators p and q the rth-degree finite-element approximation is given by 
equations (4.1). While we could determine the coefficients on the nth interval from those 
on the (n — l)st interval, attempting to determine the coefficients on the first interval, even 
in principle, leads to an inconsistency. This is because the coefficients ajt are operators. 

We illustrate this problem by a simple example for which r = 2. On the first interval 
we represent 

X 



p{x) = pq + ai- + a2j^ (4.2a) 



and 



2 

X , X 



q{x) =qo+bi- + b2-^. (4.2b) 



For the sake of complete generality we impose the differential equations (2.11) at ah and 
Ph, respectively, where < a < P < 1 are as yet undetermined: 

^1 I . I 2 /A Q \ 

— + —r-oi = Po + + O2C1; , (4.3a) 

ft lb 

^ + = -V\q, + 61/3 + 62/3'). (4.3b) 
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Next, we impose the initial conditions. The condition (2.12) reads 



[qo,Po]=i. (4.4) 

We also impose two more commutator conditions which we obtain from the continuum 
equations at t = 0. 

[Q(0),p(0)] = [Qo,ai]=0, (4.5) 
[g(0),p(0)] = [&i,po]=0. (4.6) 
There are two more commutator conditions that follow from the equations of motion (4.3): 

[6i + 2q;62,Po + ctai + «^a2] = 0, (4.7) 

[ai + 2Pa2,qo + Pbi + pH2]=0. (4.8) 

The five commutators (4.4-8) are kinematical in nature; they make no reference to the 
dynamical content of the theory, which is embodied in the function V. 

For this quantum system to be internally consistent, (4.4-6), the analogs of the three 
equal-time commutators, must hold again a,t t — h; that is 

[qo + h + b2,Po + 0'i + a2] = i, (4.9) 

[go + ^>i + ^2,ai + 2a2] = 0, (4.10a) 

[po + oi + a2,6i+262] =0. (4.10b) 

We can show that if (4.10) is assumed to hold, then (4.9) holds so long as a + P = 1. 
However, (4.10) does not hold in general unless a — 1 and (3 — 1, which implies the failure 
of (4.9). 

This kind of demonstration can be given for any stiff approximation to a quantum 
system. Thus, on the basis of quantum-mechanical consistency we reject any kind of stiff 
finite-element scheme in which more than a single initial commutator is imposed. 

Also, even if a successful stiff scheme could be found, we would prefer not to use it 
because stiff schemes are not as accurate as floppy approximations. A maximally stiff ap- 
proximation yields a relative error of N~'^ between the exact solution and the finite-element 
approximation to the exact solution. On the other hand, for a floppy approximation, the 
relative error between these two quantities is N~'^^. 

B. Consistency of the Floppy Approximation 

The failure of the stiff approximation discussed the previous subsection is not very 
surprising and it is all the more remarkable that the floppy approximation is successful. 
We begin by examining the case r = 1, using the notation of the preceding subsection. 

1. Case r = 1 

Imposing equations (2.11) at t = ah and t = Ph, respectively, yields [here qi = q{lh) 
and pi = p{lh)] 

^^—j-^ = {l-a)po + api, (4.11) 
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?l-^ = -V'[{l-f3)qo + f3q,l (4.12) 
in the first finite element. Since p(0) = po and g(0) = we have 

oo=Po, bo = (lo 
and the continuity conditions ok, t = h are 

ai^Pi-Po, ^>i=gi- go- 
Is there a choice for a and /? such that (4.11) and (4.12) together with (2.11), the equal- 
time commutator at t = 0, imply that [gi,pi] = P. Equations (4.11) and (4.12) yield the 
following commutators: 

[qi - Qo, (1 - Oi)Po + Oipl] = 0, (4.13) 

[(1 -/3)go + Pi -Po]=0. (4.14) 

Combining the three commutation relations (2.11), (4.13) and (4.14) does indeed yield 
[gi)Pi] = * provided that a and P satisfy the constraint 

a + P = l. (4.15) 

Having shown consistency with quantum mechanics on the first finite element it follows on 
all finite elements by virtue of the continuity condition on p(t) and q{t) at the boundaries 
of adjacent finite elements, t = nh. 

In this section we are primarily interested in the symmetric choice a — (3 = 1/2, where 
the equations of motion are imposed at the midpoints of the finite elements. Any other 
choice for a and P breaks time-reversal symmetry and leads to numerical approximations 
which are not as accurate as in the symmetric case. 

2. Case r = 2 

Here we impose the equations of motion (2.10) twice on each finite element. In view of 
our above remarks with regard to symmetry and numerical accuracy, we will restrict our 
attention to the symmetric case where both of the equations of motion are imposed at the 
same points x = a\h and x = a2h. 

On the n = 1 finite element we have 

o 

X X 

p{x) =po + ai- + a2j-2, 

X 

q{x) = qo + ^^"^h?' 
Imposing (2.11) at a; = aih and x = gives 

h ^2b2 ^ ^2 

— + -;— Q!i = Po + aiCKi + a2Q!i, 

h h 
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G-l 2a2 ,^// , , 2\ 

— + = -V [Qo + hai + ()2ai), 

bi 2^2 2 
— + ——a2 = P0 + 010:2 + 0202) 

h h 

^ + '^a2 = -V'{qo + bia2 + b2al). (4.16) 



From these equations we obtain the kinematical commutators: 

[po + oiCKi + 020:?, bi + 2620:1] = 0, 
[ai + 2a2Q;i, go + biai + 6201] = 0, 



[po + oi«2 + CL2al, bi + 26202] = 0, 

[ai + 2a2Q;2, go + ^10:2 + 620:2] = 0- (4-17) 



By adding these commutators we can prove that 

[gi,Pi] = [go,Po] 
if and only if oi and 02 are given by 



This is the condition for quantum consistency. 
3. Case r — 3 

Now we impose the equations of motion three times on each finite element at a; = aih, 
X = 02/i, and X = a^h. Taking 

9 3 

X X X 



p{x) =pQ + ai- + a2j^ + asj-^ 



and 



9 3 

X , X , X 



qix) = qo + bi- + b2-^ + bs^^ 



we obtain six equations analogous to (4.16) from which we derive six commutator con- 
ditions analogous to (4.17). Once again we add the six commutator conditions together. 
However, now the two commutators at a; = 02/2. are weighted by the factor 8/5. The 
condition for quantum consistency is now 

01 = ^ - v/3720, 
1 
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(4.19) 



4. General Case 

The sequence of points a at which the operator equations of motion must be imposed, 
1/2 for r = 1, 1/2 ± 1/v^ for r = 2, 1/2 and 1/2 ± ^3/20 for r = 3, fits a well-known 
pattern. These numbers are the zeros of the rth Legendre polynomial Pr {2a — 1). The 
first three such polynomials are 

Pi(2q!-1) = 2a-l, 

P2(2a-1) =6a^-6a + l, 

P3{2a - 1) = {2a - l){10a'^ - 10a + 1). 

These zeros are the so-called Gaussian knots or nodes which are used to perform quadra- 
ture integration. The weighting of the commutators necessary to derive the consistency 
condition (the factor of 8/5 mentioned above) is exactly the weighting used in Gaussian 
quadrature. -"^^ 

We conclude this section by reemphasizing that the only way to preserve the equal- 
time commutation relations is to impose the operator equations of motion at the Gaussian 
nodes. If the commutator relations are preserved at successive intervals of time, then the 
theory is unitary; that is, there exists a transfer operator [like that in (2.21) for linear 
finite elements] that is unitary and therefore probability is conserved. This same point 
has been observed in a totally different context by Durand,^^ who showed that a lattice 
discretization of the Schrodinger equation preserves orthonormality of the wave functions 
only if the lattice points lie at the Gaussian knots. 
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V. OTHER ALGEBRAS 



The successful discretization of operator equations arising from Hamiltonian systems 
associated with the Heisenberg algebra (2.12) raises an obvious question; namely, is it 
possible to discretize the operator equations associated with other algebras? For example, 
consider the algebra associated with the rotation group SO{3): 



[X, Y] = iZ, [Y, Z] = iX, [Z, X] = iY. 



(5.1) 



If we construct a Hamiltonian, H{X, F, Z), associated with this algebra, then the contin- 
uum equations of motion for the operator variables X{t), Y{t), and Z{t) are 

X{t) = -i[X,H], 



Yit) = -i[Y,H], 
Z(t) = -i[Z,H]. 



(5.2) 



The exact solution to the operator equations (5.2) satisfies the ETCR in (5.1). 

It is not easy to find a way to approximate such a system by a set of operator difference 
equations that preserve the commutators in (5.1) at each time step. For example, consider 
the Hamiltonian spin system 



H{X, Y, Z) = XY + YX. 
For this system, the continuum equations of motion are 



(5.3) 



X{t) = xz + zx, 

Y{t) = -YZ - ZY, 
Z{t) = 2Y^ - 2X^. 



(5.4) 



The linear finite-element prescription for discretizing continuum equations of motion re- 
places undifferentiated variables by averages and first derivatives by forward differences. 
Here we consider two adjacent lattice sites which we label 1 and 2. In terms of the discrete 
X, Y, and Z variables the equations of motion (5.4) become 



X2 — Xi 
h 

Y2-Yi 
h 



X2 + Xi \ ( Z2 -\- Z\ \ ^ ( Z2 -\- Z\\ ( X2 + Xi 



Z2 — Zi 

h 



= 2 



Y2 + Yi\ f Z2 + Z_ 



Y2 + Y^ 



Z2 + Zi\ fY2 + Y, 



(5.5) 
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where h is the lattice spacing. The statement of unitarity is that the SO (3) commutation 
relations hold at each lattice site. That is, if Xi, Yi, and Zi satisfy 

[Xi, Yi] = iZi, [Fi, Zi] = iXi, [Zi, Zi] = lYi, (5.6) 

then as a consequence of (5.5), X2, Y2, and Z2 satisfy the same equations. 

It is easy to see that Eqs. (5.5) violate unitarity. We merely solve (5.5) perturbatively 
in powers of the lattice spacing h. We seek solutions for X2, Y2, and Z2 in the form 

X2 = Xi + Ah + Bh^ + Ch^ + ..., 

Y2 = Yi + Dh + Eh"^ + Fh^ + . . . , 

Z2 = Zi + Gh + Hh? +Ih^ + ..., (5.7) 

We insert (5.7) into (5.5) and compare powers of h to obtain explicit solutions for the 
operator coefficients A, B, C, The results are 

A^XiZi + ZiXi, 

B = 2Y1X1Y1 + 2Z1X1Z1 - 2Xf + -Xi, 

C = Z^XiZi + ZiXiZ^ - A{XlZtXt + XiZiXl) - ^{XiZi + ZiXi), 

D = -Y,Z^ - Z,Y^, 
E = 2Z^Y^Zx + 2XiyiXi - 2Y^ + ^Yi, 

F = -ZlY^Z^ + Z^Y^Zl + 4(y2Ziyi + Y^Z^Yl) + \{YxZ^ + Z^Y^), 

G = 2{Yl-Xl), 
H = -AYiZiYi - 4X1Z1X1 - 2Zi, 

/ = 4{Xt - Y^^) + GYiZfYi - GXiZfXi + ^(y^ - Xf). (5.8) 

We can now compute the commutator of X2 with Y2 to order h^: 

[X2, Y2\ = iZ2 + Aiiy^zlY^ - XiZlXi)^ + t{xf - y2)/i3 + o{h^). (5.9) 

Observe that the difference between [X2,y2] and iZ2 is not zero but rather of order h^, 
showing that the lattice equations (5.5) violate unitarity. It is not surprising that the 
violation first occurs in order because we already know that for any differential equation 
the expansion (5.7) agrees with the continuum result through order /i^, and that the 
continuum equations, of course, do not violate unitarity. We do not know of any simple 
way to avoid this violation of unitarity. 
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Nevertheless, there is an indirect technique for obtaining a unitary set of lattice equa- 
tions of motion. We merely convert the spin system to an equivalent Heisenberg system 
having two degrees of freedom by means of the Schwinger transformation: 

^ = ^K^i - m){Q2 + ip2) + (g2 - ip2){qi + ipi)], 

Y = -^[-{Qi -iPi){Q2+ip2) + {q2 - ip2){qi + ipi)], 

^ = \[{Qi - m){Qi + m) - {Q2 - ip2){q2 + ip2)]- (5.10) 

The sorts of equations that result from this transcription will be discussed in the next 
section. 
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VI. SYSTEMS WITH MANY DEGREES OF FREEDOM 



For systems having more than one degree of freedom it is not as easy to show that the 
method of finite elements is consistent with unitarity. We iUustrate this situation with a 
system having two degrees of freedom {p, q) and (tt, (/>). For the Hamiltonian 

iy=y + Y + n5,0), (6.1) 

Hamilton's equations in the continuum read 

g = (6.2a) 

= TT, (6.2b) 
d 

p=- — V{q,cl>), (6.2c) 



^=-|r^(«'<^)- (6-2d) 



dq 
d_ 

On the lattice, the linear finite-element transcription of these equations of motion is 

(6.3a) 
(6.3b) 

qi + qo 0l + 0o\ ^g_3^^ 





- Qo 


Pi +Po 




h 


2 ' 


01 




TTi + TTo 




h 


~ 2 ' 


Pi 


-Po 






h 












h 





2 ' 2 

<?i + qo 01 + 00 



(6.3d) 



Note that while there is no operator-ordering problem in (6.2c) and (6. 2d) because 
[q{t), 0(t)] — 0, there appears to be a serious ordering problem in (6.3c) and (6.3d) because 
it is not clear that {qi + qo)/2 and (</>! -|-</)o)/2 commute. To resolve this problem we define 

2po 4go ^ 27ro 4(j)o , . 



and 



go + gi 00 + 01 

(7= — - — , T= . (6.5) 



Now, (6.3c) and (6.3d) become 



d 4 

a = —V(a, t) + -a, (6.6a) 
da h 

/5=|:^(^'^)+^^- (6-6b) 
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The simultaneous solution of (6.6) has the form 

a — a{a,(3), T — T{a,(3), (6-7) 

But, a and P involve only operators at the initial time, so [a, (3] = 0. Thus, a and r also 
commute, and there is, in fact, no ordering problem in (6.3). 

It is now necessary to verify that the ETCR's are preserved in time. From the solutions 

qi = -qo + 2a{a, P), (6.8a) 
4 4 

Pi = -Po- j^qi + j^<7{a,P), (6.8b) 

01 = -0o + 2r(a,/3), (6.8c) 
4 4 

71"! = -TTo - -01 + ^'^(«) P)^ (6.8d) 

it is necessary to verify that 

= [0i,7ri] = z (6.9a) 

and 

[qiAl] = kl^TTi] = [01, Pl] = bl^TTl] = 0. (6.9b) 

It is easy to verify (6.9a): 



4 

[qi,Pi] = [qo,Po] + -^[qo,(j{a,P)] + 2\po,a{a,P)] 

_ . 4da2i ^da / 4i 
hda h da\h^ 

= i. 



(6.10) 



However, it is harder to verify (6.9b): 

[q'i,0i] = -2[(?o,r(a,/3)] - 2[c7(q;, /3), 0o] 

AidT Aida (6.11) 
h da h dp' 

We can show that [q'i,0i] = by verifying that dr/da = da/dp because the system 
defined by (6.3) is Hamiltonian. Explicitly, we have 



dr da 



dcrdi 



da dp (0 + ^)(0 + ^) o^v 



(6.12) 



dadi 



The other commutators are evaluated similarly. For details and discussion of the time 
evolution operator see Ref. 18. 

We conclude this section with a discussion of a coupled fermion-boson system. Consider 
the system governed by the continuum Hamiltonian 

2 

H + S{x)+^'^W{x). (6.13) 
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The continuum Heisenberg equations of motion for this system are 



X = p, (6.14a) 

p ^ -S'{x) - ^W'{x)ij, (6.14b) 

= W{x)ij, (6.14c) 

-0 = ^W{x). (6.14d) 

To derive these equations wc make use of the canonical equal-time commutation and an- 
ticommutation relations for the dynamical variables 

[x, p] = i, (6.15a) 

[V'>]+ = 1, _ _ (6.15b) 

[x, t/j] = [x, ij] = = [p, -0] = 0, (6.15c) 

^2 = ^2 ^ 0. (6.15d) 



In (6.14) we have ordered the operators in order to make the equations of motion manifestly 
Hermitian, which ordering is irrelevant in the continuum, but not on the lattice. 

If we were now to put the system (6.14) on the lattice using our finite-element pre- 
scription, we would find that the resulting lattice theory is not unitary. In anticipation of 
this difficulty we will replace the function W' in (6.14b) by a function Z, which will be 
determined by the requirement that the system be unitary. We will see that Z satisfies an 
interesting nonlinear equation. We will, of course, find that as the lattice spacing h tends 
to zero, Z approaches W'. 

We take, then, for our lattice difference equation 

(6.16a) 



Xi 


- Xo 




h 


Pi 


-Po 




h 




- i'o 




h 




--00 


h 



2 ' 

-S'{a) - ^Z{a)(i), (6.16b) 

W(a)(f), (6.16c) 

'^W{a), (6.16d) 



where 



a.^^, (6.17a) 

(6.17b) 

We will now simply quote the results found in Ref. 18. There is no problem with 
unitarity in the pure fermion sector of the theory, because it is easy to show that 

l-ihW(a)/2 , 
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Similarly, there is no problem in the pure boson sector. The constraint comes when we 
examine the mixed commutator [a;i,7/'i]. Requiring this to vanish leads to the following 
equation for Z: 

Zix) = 



+ [1 + h^W\x)/4] (^S'{x) -S' (^x--^ 



h^W'{x)/A 



+ h?W'^{x)/A){l + h?S"{x)/A) 

(6.19) 

For small lattice spacing this approaches 

zw.w'(x)-^::m«^+o(/.«). (6.20) 

Notice that if /S is a polynomial of degree two, Z is exactly W . This suggests that 
electrodynamics will not present subtleties when analyzed using finite elements, while 
super symmetry, for example, may require more care. 
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VII. SCALAR QUANTUM FIELD THEORY 



It is straightforward to generalize the above discussion to quantum field theory. In this 
section we will apply the method of finite elements to self-interacting, two-dimensional 
scalar field theory, and in particular, calculate the mass renormalization for the {(jP'^)2 
and the sine- Gordon field theories. 

Consider the Hamiltonian 

H=\<i>U\<i>l + V{ct>). (7.1) 
This gives rise to the operator Klein-Gordon equation 

ht-^xx = f{4>), (7.2a) 

where 

/(0) = -V'{4>). (7.2b) 

Because we will be using linear finite elements, we rewrite (7.2a) as a system of first-order 
equations: 

(j)t = TT, (7.3a) 
4>x = r, (7.3b) 

7Tt-r, = f {<!>). (7.3c) 

We introduce a rectangular finite-element lattice with the time-lattice spacing being h and 
the space-lattice spacing being A. The spatial extent of the lattice is L. If we approximate 
the field in the finite element by a polynomial linear in x and t, 

(l){x,t) = (^1-^ (l- (f)m-i,n-l + ~ ^0m,n-l 



and impose the equations of motion (7.3) at the center of the finite element, we obtain the 
system of difference equations 

^(0m,n+l + 0m+l,n+l — <t>m,n — <t>m+l,n) 

= ^ i'^m+l,n+l + TTm.n+l + '^m+l,n + '^m,n), (7.5a) 
2^(0m-|-l,n-|-l + 4>m+l,n — 0m,n-|-l — 4'm,n) 

+ + J- m+l,n 

), (7.5b) 

'^J^{j^m,n-\-'l ~l~ TTm -I- l,n -1-1 '^m,n '^m+l,n) 0-^m+l,n+l ~l~ ^m+l,n ^m,n+l ^m,n) 

— f ^rn+l,n+l + 0m,n-|-l + 4>m+l,n + 0m,n j ^^-^ 
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The first index, m, on the fields represents the spatial lattice site and the second index, 
n represents the temporal lattice site. Note the significant mnemonic for linear finite 
elements: corresponding to directions in which derivatives are taken, a forward difference 
is taken, while in other directions, a forward average is taken. The index m ranges from 
to M = L/A; all fields are taken to be periodic in the spatial lattice: 



50,n = (PM,n- 



(7.6) 



For technical reasons, which will be discussed below, we will take M to be odd. As shown 
in Ref. 19, it is the space-averaged operators, 

^m,n = -j^{(i>m,n + 0m-l,n), (7.7a) 
nm,n = 2 + 7rm-l,n), ('''•7b) 

which obey the canonical equal-time commutation relations 

[^m,n)^m',n] — [-n}7i,n) n^'^^] — 0, 

These commutation relations are the discrete analogs of the continuum equal-time com- 
mutation relations. 

In principle, we can solve the system of difference equations (7.5); for the purposes 
here, however, it is sufficient to expand in powers of the temporal lattice spacing h. We 
do not expand in A. The expansions for (f)m,n+ij '^m,n+ij ^ind F^^^+i are 

4>m,n+l = 4>m,n + hAm + h^Bm + 0{h^), 

+ hCm + h^Dm + 0{h'^), (7.9) 

Inserting (7.9) into (7.5) leads to a set of difference equations for the operators A^, B^, 
Cm, These equations all have the same generic form: 

Xm+Xm+l=Rm, (7.10) 

which, for periodic boundary conditions and M odd, has the general solution 
1 



M-l m-1 



M+m-1 

.k=m fc=0 



1 J2 {-If+^R,. (7.11) 



2 

k=m 



Next, following the quantum-mechanical discussion in Sec. II, we introduce a Fock-space 
representation for the canonical operators $rn,n ^ind 11^ „, deffned in (7.7): 



M 

fc=l 

M 

n-,n = J] + aie-^'^^^VM)^ (7_;L2b) 
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where 

[afc,aj] = 5k,h (7.13a) 
[ak,ai] = [alal]=0. (7.13b) 

In (7.12) 7fc, k = 1,2,...,M, are arbitrary parameters which later wiU be fixed by a 
variational argument, similar to that given in Sec. II. Recall that the spatial size of the 
lattice is L = MA. 

As a consequence of (7.13), the equal-time commutation relations (7.8) are satisfied at 
the initial time n provided that 

lk = lh-k- (7.14). 

It is crucial that as a consequence of the operator difference equations (7.5) the equal time 
commutation relations (7.8) are preserved at all subsequent lattice sites, in particular those 
at n + 1. 

Using (7.12) we find easily 

2(im-i,n + im,nj - -^2^7/= tan— (OfeC ' - a^e ), (7.15) 

fe=l 

It is easy to verify the expected equal-time commutation relation between ^{Tm-i,n+^m,n) 
and ^m,n- 

m— l,n + Tm,n)] = 0. (7.16) 

We will merely display the order- /i coefficients: 

Am = lim, (7.17a) 

UEm-i + Em) = ^Y.- tan ^(afee^'^-^-Z- + ale-'"^'^/^), (7.17b) 

fc=i ^'^ 

1 A ^ k- 

^(Cm-i + Cm) = -^Y.^kt^''" |^(a.e^'^-2-/- + aie-^'^^'-n 

k=l 

+ fi^m,n). (7.17c) 

We do not bother to display the solutions for Bm, Dmj and F^; however, they can be 
calculated similarly. 

Consider first a free theory, for which / = — /i^^. We extract spectral information by 
taking matrix elements of (7.9) between the Fock vacuum and a one-particle state with 
lattice momentum I: 

(l,^| = (0|a/. (7.18) 

We compare these matrix elements with an assumed approximate exponential time depen- 
dence with a single frequency: 

{l,l\^m,n+l\0) ~ e^"'^l,/|$m,n|0), 

1 (7.19) 
{l,l\Um,n+l\0) - e"-^''{l,l\Um,n\0)- 
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The 0{h^) and 0{h^) equations give a relation between the frequency uji and the variational 
parameter 7/, 

(7.20) 



LOl = 



and the dispersion relation 



2 2,42^^ 



(7.21) 



In the continuum limit, A — > 0, M — > 00, and L = MA — > 00, the lattice equivalent of the 
continuum momentum p is 



P = 



2ttI 
MA" 



Thus, we recover the continuum dispersion relation 

2 2,2 

UJ = jJL +p . 



(7.22) 



(7.23) 



It is remarkable that if we include the two additional equations coming from the 0{h^) 
terms in $^ ^+1 and Ilm,n+i as well as the three further equations coming from r^^^+i 
only redundant information is supplied. 

Consider now an interacting theory for which 



2 2N^ 



(7.24) 



Following the procedure described above, we find that (7.20) still holds, but the dispersion 
relation (7.21) is replaced by 



2 2,4 2^^ 

=TOren + ^tan — , 



where 



M 



A.^ + (2iV-l)!!^,X^-\ X = Y,^l 



(7.25) 



(7.26) 



fe=i 



In the continuum limit, an asymptotic analysis^° serves to evaluate X 



X = — In , 

27r V m™ A 



(7.27) 



An equation for the renormalized mass is thus obtained when we substitute (7.27) into 
(7.26): 



2 , (2iv-i)n^ 



(27r) 



N-l 



In 



m™A 



N-l 



(7.28) 



This nonperturbative result closely resembles the formula in continuum perturbation the- 
ory: 



rrir. 



{2N -l)\\g 
(27r)^-i 
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In 



n AT-l 



(7.29) 



The correspondence is provided by the identification of the momentum cutoff A with tt/A. 
We close this section by extending this calculation to the sine-Gordon model, for which 



V =^cos gcf). (7.30) 
9 



The corresponding force is 



2 °° (_^\N 2Nj^2N+1 

f = sin# = -fx' y ^ V . • (7-31) 

Each term in the sum in (7.31) gives a contribution to the renormalized mass of the form 
given in (7.26), so the formula for the renormalized mass is therefore 



2 2 V- (-1) r^(2A^ + l)!!^ 2 

Ar=o ^ 

where X is given in (7.26). The dispersion relation is 

nl 



An asymptotic analysis in the continuum limit leads once again to (7.27), so a simple 
calculation yields the following relation between the renormalized and the unrenormalized 
masses: 

TOren = A* ( — I • ('•34) 

This is the characteristic power-law renormalization found in the conventional treatments 
of the sine-Gordon model^^, and reduces to the perturbative result of Coleman^^ when 
g'^/Sir is small. 
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VIII. THE DIRAC EQUATION AND FERMION DOUBLING 

The finite-element lattice Dirac equation in 3 + 1 dimensions is 

-^(^m,n+l - 1pm,n) + ^ (V'm^ + l,m_L ,n " V^m,-,m_L,n) + /^V'm.n = 0. (8.1) 

Here the overbar represents a forward average over that coordinate: 

__ 1 . 

and the notation TTix means that all spatial coordinates but rrij are averaged. Let us begin 
by finding the momentum-space spinors, the eigenvectors of the transfer matrix. That is, 
write for a plane wave at time n 

and at time n+1 

Vfe,„+i = Un+ie-'^-"'^^'''. (8.2b) 
The transfer matrix T is defined by 

ttn+l = TUn- (8.3) 

By substituting (8.2a) and (8.2b) into the Dirac equation (8.1) we easily find that 

■ - + \ -1 

!2_ + 111 + 

h A 2 



h 


7-t 
A 




(- 


4 


A2 



,1+^ + ^'') (^l-^-^<^ + in"7-t + MW7''j. (8.4) 
where 

t = tp, (tp)i = tanpj7r/M. (8.5) 
Let us adopt a representation of the Dirac matrices in which 

^°={l -1)' = :). (8.6) 



Then the eigenvalues of T are easily found: 

1 

1 ^ ihui/2' 

Here Co is an abbreviation for 



l±z/ic2;/2 , , 

^=V—T-^- (8-7) 



'4^2 
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which is exactly the same as the dispersion relation (7.21). It is obvious that A has modulus 

unity, so it can be written in the form A = exp{±iu)h), where u is, of course, a function of 
h. (The relation between uj and u) is = tan ^.) The corresponding eigenvectors may 
also be found straightforwardly. They are to be normalized according to 

w±7°w± = ±1. (8.9) 



They are 

where x is a two-component, rest-frame spinor, normalized by y^x — 1- Thus, with 

n5=(j J), (8.11) 



we have 

\ / ~ . \l/2 ^,/_ X 1/21 

,(0) 



~ , \ 1/2 - + / ~ \ 1/2' 

±^75 



<^ (8.12) 



where is a four-component rest-frame spinor with 7*^ eigenvalue of ±1. Therefore, in 
terms of the spinors «+(?) = ^^-(p), W-(p) = u_(— p), we have the completeness relations 

«±wi7° = ±^ (^// ± 7°^^ T ^) , (8.13) 



spins 



which in the continuum limit reduces to ±(/i ^ 7 • p)/2//. We have the same result on the 
lattice, provided we define 

2t 

p=-. (8.14) 
All of this tells us that the momentum expansion of the Dirac field has the form 

= E \/| (^v^su^^s^^-^-"" + dl,v^,,e-'^-^-l'') , (8.15) 

where we now use the standard notation u = i^^U-, v = i'j^u^, with the usual inter- 
pretation that cP creates a positive energy positron, while b annihilates a positive energy 
electron. The canonical lattice anticommutation relations 

{'ipm,n,'iplii'^J ^ (8.16) 

will now be satisfied if 

{^P,s'^p',s'} = j3^P,p'^s,s', {dp,s,dl,^^,} = -^(5p,p'(5,,,/, (8.17) 
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and all other anticommutators of these operators vanish. 

The unitarity of T is sufficient to establish the unitarity of the fermion sector in the 
noninteracting theory. For further details see Refs. 23 and 24. 

It is apparent that the dispersion relation (8.8) solves the fermion doubling problem; 
that is, oj^ assumes the value /x^ at p = (mod M) and nowhere else.^^ In this discussion 
we have assumed, as in Sec. VII, that M is odd so that the Dirac field is periodic; however, 
the same conclusion would hold if antiperiodic boundary conditions were used. In the 
dispersion relation one would simply replace phy p + ^. 

Let us conclude this section by summarizing the properties of the linear finite-element 
Dirac equation (8.1): 

1. It is unitary in that the equal-time anticommutation relations are exactly preserved in 
time. 

2. It may be derived from an Hermitian action. 

3. There is no fermion doubling. 

4. The difference equation is local, in that only nearest-neighbor terms appear. 

5. It is chirally symmetric in the massless limit. 
For a complete discussion of these points see Ref. 23. 

The no-go theorems of Karsten and Smit, Nielsen and Ninomiya, and Rabin^^ are 
avoided because the time development operator is nonlocal, which arises because undiffer- 
entiated fields appear as averages. 
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IX. DIRAC EQUATION WITH INTERACTIONS AND THE SCHWINGER 
MODEL 

In the last two sections of this review, we wiU discuss interactions of fermions with 
gauge fields. In this section we will concentrate on the Dirac equation, interacting with 
either Abelian or nonabelian gauge fields, and make application to electrodynamics, in 
particular to the Schwinger model. In the next section the nonabelian interactions of the 
gauge fields among themselves will be derived. (The original treatment of finite-element 
electrodynamics was given in Ref. 27.) For simplicity, initially our discussion will be 
restricted to (1+1) dimensions. 

A. Equations of Motion in the Continuum 

We begin by recalling the form of the continuum field equations of a nonabelian gauge 
field coupled to a fermion field ip. Let us start with the free Dirac equation 

{0 + 1^)^ = 0. (9.1) 

Equation (9.1) is invariant under an infinitesimal gauge transformation 

i/j ^ i/j + Si/j, Si/j = igSuji/j, Suj = SujaTa, (9-2) 

provided 6u> is constant. Here T is the generator of the gauge group. If du is not constant, 
we can restore the invariance by adding an interaction term to the Dirac equation, 

i'i$ + g^ + fi)il; = 0, A^' = A>'^Ta, (9.3) 

provided A transforms according to 

A^^ A^ + 5A^, 5A^ = d^Suj + ig[Su, A^]. (9.4) 

Under (9.4) the curl of A is not covariant, so we must add a suitable interaction term to 
construct the field strength: 

i^M^ = di,A^ - d^Af, - ig[Af,, A^]. (9.5) 
The field strength transforms covariantly: 

+ SF/^jy, 6F^jy = ig[Su;, F^^]. (9.6) 

Finally, because the current 

r = g^T^^^ (9.7) 

transforms covariantly, 

r^r + Sr, 5r = ig[5u;,j% (9.8) 
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we must include the interaction term in the Yang-Mills equation, 

D^Fi''' = f, where = - ig[A^, ]. (9.9) 

B. Free Lattice Equations of Motion 

In 1 + 1 dimensions, the free Dirac equation (8.1) reads 

-J^{4'm,n+l — 4'm,n) + -^(^m+l,n — ^m,n) + -^{Om+l,n + ^m,n) = 0. (9.10) 

Here we have used the abbreviations 

<i)m,n = ^(V'm+l.n + '0m,n), (9.11a) 
Om,n = ^(V'm.n+l + V'm.n)- (9.11b) 

Similarly, the field strength E = Fqi is constructed as 

-^m,n = ^(Cm+l,n+l " C'm+l.n) - -^{Bm+l,n+l " -Bm,n+l), (9-12) 

where 

Bm,n = 2[(^o)m,n + (^o)m,n-l]5 (9.13a) 

Cm,n = ^[(^l)m,n + {Ai)m-l,n]- (9.13b) 

In (9.12) the tilde signifies the average over the four adjacent lattice sites: 

and the (0) superscript is a reminder that this is the free field strength. 
Finally, the free Yang-Mills equations driven by a current are 

l{F^+l,n - i^,n,n) = "JI,., (9.15a) 

-^{Gm,n+1 - Gm,n) = Jm,n' (9.15b) 

where 

Fm,n = ■^iEm,n+l + Em,n), (9.16a) 

Gm,n = 2^Em+l,n + Em,n), (9.16b) 
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and 

3rn,n ~ 9^ m,nT'y^^ rn,n- (9--'-'^) 

Here we have used = ip. (We have, of course, anticipated the interaction with the 
fermion by including j^.) 

C. Constructing the Interaction Term Involving 

We now proceed to construct the interaction terms for the Dirac equation (9.10). We 
begin by recognizing that the gauge transformation of the fermion on the lattice is 

(^*m,n = igSu!m,n^m,n, (9.18) 

which is the appropriate lattice version of (9.2) because (9.18) guarantees that the mass 
term in (9.10) transforms covariantly, that is, by the same rule. (We also note that then 
the current in (9.17) transforms covariantly.) We can regard (9.18) as a difference equation 
for 5'^rn,n'i in particular, we find 

M+m 

50m,n = ig5uJm,n0m,n + ^ XI (-1)™+™ {SuJm' ,n " 5uJm' -l,n)Om' ,n- (9.19) 

m'=m+l 

Here, we have used either periodic or antiperiodic boundary conditions depending on the 
size of the lattice: 

V^m+M,n - (-1)^+Vm,n, #m+M,n = (-l)^+'5Vm,n, (9.20) 

where M is the number of spatial lattice sites. The boson variables such as 5ujrn,n will be 
assumed to be periodic. 

It is evident that the second term in (9.10) does not transform covariantly under (9.19); 
as in the continuum, an interaction term must be added. Now the free field strength in 
(9.12) is invariant under the transformations 

{5AQ)ra,n = ^(^Am+l,n+l + 5Am,n+l - ^Am+l,n " (^Am,n), (9.21a) 

{^M)m,n = ^^{^^m+l,n+l + ^^m+l,n " 5Kra,n+l " 5Km,n), (9.21b) 

The finite-element connection between 5Am,n and 5ujm,n is 

m+l,n+l 

Then the scalar and vector potentials, (9.13a) and (9.13b), transform by 

(^^°^-Bm,n = J^{Su;m,n " ^'^m.n-l), (9.23a) 
S^^^Cm,n = ^{S0Jm,n - (^"^m-l.n)- (9.23b) 
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Here the superscript (0) reminds us that further transformations of B and C will have to 
be deduced. 

In the rest of this subsection wc will examine that portion of the Dirac equation pro- 
portional to 7^. A short calculation reveals that under (9.19) 



^ (^m+l,n ^m,n) 



m+M 



E (-1) 



m+m' 



S^m'-l,n)(^m',r 



m'=m+l 



. (9.24) 



The first term here expresses the desired covariance of this term in the Dirac equation. 
The second term will be cancelled if we introduce an interaction term 



m+M 



(9.25) 



m'=m+l 



and vary it with respect to C according to (9.23b). Indeed, in the continuum limit, this 
term reduces to the corresponding interaction term in (9.3) because 



m+M 

E (-1) 

m'=m+l 



m+m 



— ^(mA) as A — > 0, m — > oo, M — > oo. 



(2.26) 



if QM+m = (— l)^^^^m- But, On the lattice, we are not finished, for we have neither yet 
varied (9.25) with respect to 6 according to (9.19), nor achieved the covariance of the 
interaction term. 

A straightforward calculation now reveals that 



"OJ-m,n ''y"^m,n-i-m,n " -^m,n " -'m,n- 



(9.27) 



Here the first term is the required covariance term, the second term involves a new variation 
of C, 

S''^^Cm,n = ig^[SuJm,n + Su)m-l,n, C'm,n], (9.28) 



and the third term is the 5^^^ variation of a new interaction term 



. I m+M m'+M-l 

m'=m+l m"=m'+l 



m ,n^m ,n- 



(9.29) 



It is easy to see that (9.29) vanishes in the continuum limit, while the variation (9.28) 
reduces to the second term in (9.4) in that same limit. Once again, we are not finished: we 
still need to vary /^^^ with respect to 56 and S^^^C, and produce the required covariance 

0f/(2). 
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Clearly, this process of adding successive interaction terms and evermore C variations 

never terminates. But it is easy to discern the general pattern. The easiest way to express 
where is the order in gAC, is by the following inductive formula: (here and in 
the rest of this section we delete the time index n, since all variables are evaluated at that 
time) 

m+M N 

N>1: 4^^ = ^ E (-ir+'"TrfH^AC^')'4^-'\ (9.30) 



m — m+l k=l 



where we define I^n} = —2i^^9m/ A- The gauge transformations are given by (9.19) and 

kj^l: S^'^^Cm = l!^^ [. . . [dcv^ - 5um-i. Cml-.., (9.31a) 

k nested commutators 

<5^'^ = f \^^m + ^^m-U Cm] , (9.31b) 

where Bfc is the A;th Bernoulli number. The required covariance statement 

AT 

hl^r^^^ + E <5(^)/(f = igSUrnlL""^ (9.32) 

is proved in Ref. 28. 

From (9.30) we can derive an "integral equation" satisfied by the full interaction term 
for the vector potential, / = Y1'n=i ^^^^ '• 

m+M 

Im = l Yl (-ir+"^'(e--^^-'-l)(4°.^+/m'). (9.33) 

m'=m+l 

From this, a difference equation can be immediately derived: 

Im + = ^ (1 - e'"^^-) (9.34) 



D. Scalar Potential Part of Dirac Equation 

The procedure is now clear. We start from the part of the free Dirac equation (9.10) 
proportional to 7^, 

-j^{4>n+l-4>n). (9.35) 

Here we have dropped the spatial index m because throughout this section all variables 
will be evaluated at the same spatial coordinate. Now, we need to solve the fermion 
transformation equation (9.18) for 50. We immediately find 

S4)n = ig5(XJn(j)n + {-'^TiHo " WSu!o(f)o) 

n 

-igJ2 (-l)"+"'(5'^n' - du;n'-i)<Pn'. (9.36) 

n' = l 
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At this point we make a slight variation on the procedure of Ref. 27. We will choose as an 
initial condition on the variation 



Ho = -^{Sujo + (5a;_i)0o, 



(9.37) 



which will simplify the form of subsequent formulas. The reason for the choice (9.37) (or 
the slightly different choice made in Ref. 27) is to ensure that the lattice variation of the 
time difference of and hence the Dirac equation have the correct continuum limit. The 
latter is obtained from 

n 

2 J](-l)"+"7„, + (-l)"/o^/(n/i), as h ^ 0, n^oo. (9.38) 

n'=l 

Using (9.37) we write (9.36) as 

n 

5(j)n = ig5uJn(j)n - (^SuJn' - 6uJn' -l)4>n' , (9.39) 

n=0 

where the prime on the sum signifies that the 0th term is counted with half weight: 

n n 

J]'/n' = ^/n' + ^/0. (9.40) 



n'=0 



n=l 



The variation of (9.35) is immediate: 



— (0n+l - 0n) 



^-^ig5Un{(j)n+l - (l>n) 



■ 

+ 2z^^ Y^i-ir+'^'iSujn' - 5u;n'-l)(f>n'. (9.41) 



n'=0 



The first term here is the required covariance of (9.35), while the second term on the 
right-hand side of (9.41) is cancelled by the variation of the following interaction term, 



■ 



n'=0 



under (9.23a), 



(9.42) 



(9.43) 



Of course, using (9.38), we find that K^^"* reduces to the appropriate interaction term in 
the continuum Dirac equation. 
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However, here again we are not finished. We must vary K^^^ with respect to S9. Doing 
so necessitates the introduction of a new field variation, 

n > 1 : 5^^^Bn = '-^[SiOn + 5un-i,Br,], (9.44a) 

n = 0: S^^^Bo=''-^[Sujo + Suj-i,Bo] + '^[6u;o-Suj-i,Bo], (9.44b) 

which reduces to the appropriate portion of the continuum variation (9.4) in the continuum 
limit, and a new interaction term, 

Q n n' 



kP = ~-^2{ighf ^\-ir+-"Br,,Bn"<Pn", (9.45) 

n'=On"=0 

which vanishes in the continuum limit. Here 

^' f> 1 " ^ 1 

l' ■ ^ fn" = -^fo + ^ fn" + -^fn'i 



n'-l 

n"=0 ~ n"=l 





fn"=0. (9.46) 



n"=0 



Once again, this iterative procedure continues indefinitely. We may again write the 
general result in terms of an inductive formula. The order N interaction term is 



iff = - E (-1)"""' E 



n'=l fe=l 

+ l(i9hr^{-irB^4'\ (9.47) 

where Kn^ — — 2z7°0^//t. In Ref. 28 we show that these transform according to the 
required law, 

N 

S,Ki^^ + E S^'^Ki^-'^'^ = igSo^nKi'^K (9.48) 

fe=0 

where S(f)n is given by (9.39) and 

ky^l: d^'^^Bn = M)!^ [. . . [Sujr. - S,], . . . , n ^ 0, (9.49a) 

k nested commutators 

S^'-^Bn = '-jiSuJn + SuJn-1, Bn], U 0, (9.49b) 



Ml: S('^Bo ^ M)!(-l)^^ [. . . [5uo - 6u i?o], . . . , B,] 



(9.49c) 



k nested commutators 



Y2 

5^^^ Bo = '-^[Sojo + Soj-u Bo] - ig^i^ujo - 5u;_i, Bo], (9.49d) 
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where, again, Bk is the kth BernouUi number. 

As before, from (9.47) an "integral equation" can be immediately obtained for the full 
interaction term with the scalar potential, Kn = Y1n=i ^n^^'- 



It 

n' = l 

+ (_l)n [e'9hBo/2_^^^j^{0)^ 

which is equivalent to the difference equation 



The full lattice Dirac equation is given by (9.10), (9.30), and (9.47): 



(9.50) 



(9.51) 



-^(^m,n+l— ^m,n) + -^(^m+l,n — ^m,n) + ^(</'m,n+l+</'m,n) + (-fm,n + -f^rn,n) — 0. (9.52) 

Note that (9.52) gives V'm,n+i in terms of fields at time n and earlier, so that this difference 
equation may be solved by time stepping through the lattice. 

E. Unitarity of the Dirac Equation in the Temporal Gauge 

For the case of Abelian electrodynamics, the Dirac equation may be written explicitly. 
(For a complete discussion, see Ref. 27.) The generalization to four dimensions is imme- 
diate, and the result may be expressed in terms of the spatially averaged electron field as 
follows:^^ 



j\Wrn,n+l ~ Vm,n) ^ 



M 



A^O^ 'T / \ — ^ (i) 



(9.53) 



where a sum over the repeated index j is understood. Here, we have adopted a temporal 
gauge, = 0, and expressed the interaction in terms of (only the jth index is explicit 
and the spatial coordinates refer to the jth direction) 



M 



2 

X exp 



m" = l 



M 



sgn(m' 



m)sgn(m'" — m")sgn(m" — m)Cn 



m"'=l 



(9.54) 
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We have used the abbreviations 

^-^■ = V^^' C = EC...' (9.55) 

and 

= (9.56) 
We can now carry out the sum over m" in (9.54): 



M 



-1 + COS j sgn(m" — m)sgn(m" — 'm')(rn" I secC 

\m" = l / 

(-1)""+"^' sin ( J2 sgn(m" - m)sgn(m" - m')Cm") sec C, (9.57) 

where 



M 



^m" = l 



1, m' > m, 

em',m = 0, m' = m, (9.58) 
—1, m' < m. 

It is obvious that a^^^ is Hermitian. 

Let us write the Dirac equation (9.53) in the form 

U^Pn+l = Vtl^n- (9.59) 

It is apparent that V = 2 — U^ so the transfer matrix is 

T = 2U-^-l. (9.60) 
The condition that T is unitary translates into the foUowing condition on U : 

U + U^ = 2. (9.61) 

From (9.53) the matrix U is explicitly 

h I 

^^m,m' ^ 5^^, m^- (9.62) 

Therefore, the unitarity condition (9.61) is equivalent to the condition that a^^^ be Her- 
mitian: 

a^^^^ = a^^\ (9.63) 

which is satisfied as noted above. 
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F. Solution of the Abelian Equations — The Schwinger Model 

Here we wish to illustrate how the equations of motion are solved, and spectral in- 
formation extracted, for the special case of the Schwinger model, electrodynamics in two 
space-time dimensions with zero bare fermion mass. We can work in the gauge Aq = 0, 
and choose a square lattice, A — h, in which case the difference equation (9.34) can be 
solved by 

</'Si+i =e*^"^--0S^i,n, (9.64a) 
<f>^^\,n+i =e-^^'^^--0fel (9.64b) 

where the superscripts denote the chirality, the eigenvalues of ^75 = 7*^7^. We take (/)m)n, 
to be the canonical fermion variables, which for free fields have the Fock-space expansion 

M 

(f)^^^ = J2 e'P^'^^{e-'P^''^v^^^a^^^ + e'^^'^^v^-^a^'^), (9.65) 
fc=i 

where for M even, pk — {2k + l)7r/(M/i), v^"^^ are the eigenvectors of 75, and 



Wlt\^i^\ = ^S,,,, (9.66) 



the other anticommutators being zero. The physical interpretation of a^i^\ a^j^^^ as cre- 
ation and annihilation operators is as follows: 

M 



for 0<A;<y-l, 4+) |0) = 4 ^ V) = 0, 



for -^<k<0, 4+)t|o)=4-)|o) = 0; 



(9.67) 



this construction then implies the correct lattice fermion Green's function. 

The only physical particle in the Schwinger model is a boson, which we denote by B, of 
mass /X. We can obtain^^ the dispersion relation for this particle in a manner analogous to 
that employed in Sec. VII by taking matrix elements of the equations of motion between 
the vacuum and a B state of momentum qi = 27rl/Mh. Using (9.15) we find 

{B,1\J(+1\0)=^{B,1\J^^^ + J^J0) 

= ^{-(e^'^''^ + l){e-"i^^ - 1) + (e^'^''^ - l)(e-^«''^ + 1)}(S, l\Em,n\0) 

^^iqi+uJl){B,l\E^,n\0), (9.68) 

as /i — > 0. On the other hand, we evaluate the current matrix element using the solution 
of the Dirac equation (9.64) and the Fock-space expansion at the initial time n = 0. That 
is, since 

0Sn+i = e^''S"-^-"+-0^;-)^_,^„, (9.69) 
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and 



we have 



M 

\ ^ ^ipkmh „{+) 



a 



k ' 



(9.70) 



M 

Q(+)tgi(Pfe/-2?fe)(m-n)/iM _|_ gi(pfe-pj./)/i 



k,k' = l 



m-\-l ,n 



_|_ ^iPkh^-ieh y^^^_^{Cm-n. + r,'r-Cm-n + l + r,r) ^-i 
_|_ Q-iPk'h^ieh'^'^^^{Crn-n+r,r-Crn-ri+l+r,r) ^iehCrn+l,n-^^(+) _ (9-71) 



We now assume that the B states are not created by fermion operators, that is, for all 
k>0, 

(5,Z|4+)'^ = 0, (9.72) 

and that the commutator of a^j^^ and Cm,n is negligible as h ^ 0. Then, using the canonical 
relation (9.66) together with the vacuum definition (9.67), we find that the matrix element 
of (9.71) is 

e 1 



k=-M/2 



X (S, /|Cj„_|_i^^ + 'Y^{Cm-n+r,r ~ C'm-n+l+r,r) |0) , (9.73) 



r=0 



where an expansion in h has been carried out. The sum on k in (9.73) is immediately 
evaluated as 2i/ sin(7r/M), while the remaining matrix element is 



{B,l\Cm,n\^) 



+ iqihy^^< 



,ih{qi-u>i)r 



r=0 



i^i - qi 



{BJ\Cm,n\0), (9.74) 



where in the last summation on r we have deleted a rapidly oscillating term ~ e ^^^^ Qi)nh 
{nh oo). Finally, from (9.12) we learn, as /i ^ 0, that 



iuJi{B,l\Cm,n\0) = {B,l\Em,n\0). 



(9.75) 



When we put (9.68), (9.73), (9.74), and (9.75) together, we obtain the desired dispersion 
relation. 



{B,l\Em,n\0) =-^—{B,l\J^X\^) 

1 uJi + qi 

2 1 ze^ 2i 



i oji + qi AM sin7r/M oji — qi 



(S,Z|C^,„|0) 



1 

2 ~3 7777r:r7i7(^'^l-^m.'^|0)' 



cuf - q^ Msin7r/M 
45 



(9.76) 



or 

,2 2 I ,.2 



u;t = qt + fi\ (9.77) 

where the mass fj, is 



Using the solution (9.64) we may also calculate the divergences of the vector and axial 
currents. It is easily seen from (9.64) that the lattice divergence of the vector current 

Jl^,n = e^l^yr^m,n (9.79) 



2 



IS 



{"d^J,nm,n = ^Sin (^^m,n) e-'^(^-+i-+i+^-"+^)/20i^^,^^^,0^,,+l + /i.C, (9.80) 

which has vanishing vacuum expectation value because 

On the other hand, this same result shows that the axial- vector current 

J^m,n = e*1,,„7V75*m,n (9.82) 

yields the axial anomaly 

iC'd^JtUn) = -e^ (^Msin^y' Em,n. (9.83) 

The same anomaly emerges as in (9.78). The details of this latter calculation can be 
found in Ref. 27. Note that the lattice anomaly e^/(Msin(7r/M)) differs from the con- 
tinuum value c^/tt by a term of order an error typical of the linear finite-element 
approximation. 
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X. LATTICE YANG-MILLS EQUATIONS 

We now continue the discussion where we left off in Sec. IX. C. For simphcity, we 
continue to work in 1 + 1 dimensions, and we wiU here assume a square lattice, /i = A. 

A. Construction of the Field Strength 

Now that the gauge transformations for C and B are completely determined, it should 
be straightforward to work out the interaction terms in the field strength, E = Fqi. We 
begin by recalling that the zeroth order construction (9.12) 

-^m,n = -^{Cm+l,n+l - Cm+l,n) - -^{Bm+l,n+l - Bm,n+l), (10-1) 

is invariant under (9.23), 

S^^^Cm,n = J^iSuJm,n " SuJm-l,n), (10.2a) 
5(°^S^,n = l{SiOm,n ' ^O^m.n-l)- (10.2b) 

However, under the S^^"* transformations (9.28) and (9.44) the variation of E^^^ is 

) C'm+l,n+l] + [Bm,n+1, C'm+l,n+l] + [Bm+l,n+l, Cm+l,n] 
— [Bm,n+1, Cm+l,n] + [C'm+l,n, C'm+l,n+l] — [Bm,n+1, Bm+l,n+l]} ■ (10.3) 

The total variation in (10.3) allows us to identify the first-order interaction term. Because 
it is already clear that the interaction terms in the field strength will be local, involving 
fields at the four corners of the finite element, let us simplify the following by introducing 
the abbreviations 

Bi = Bm+i,n+ii and Bq = Bm,n+i, (10.4a) 
Ci = Cm+i,n+i, and Cq = Cm+l,n- (10.4b) 

Then, the free term (10.1) is 

^ffn = ^(Ci - Co - Si + 5o), (10.5) 
and the interaction deduced from (10.3) is 

^ffn = -{{[Bi, Ci] + [Bo, Ci] + [Si, Co] - [Bo, Co] + [Co, C,] - [Bo, Si]}. (10.6) 

Note that E is antisymmetric under the interchange S <-> C, an antisymmetry that will 
be maintained in each order. We could add here an arbitrary multiple of (£^^,n)^ because 
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this is invariant under 5^^^ variations. We choose not to do so, and thereby keep a minimal 
form in terms of nested commutators. This form, of course, guarantees hermiticity. 

It is straightforward to find the higher-order interaction terms. We find E^"^^ from the 
equation 

S('^E(^l + 5^'^&^l + 5('^E(^l ^ ^g[Su;m,n, ^ffj- (10.7) 
Up to the above-mentioned ambiguity it is 

^Sn = -^{[^1' [^I'^i]] - [C*!' [Ci,B^]] - [Si, [Si, Co]] + [Ci, [Ci,So]] 
-f 2[Si, [So, Ci]] - 2[Ci, [Co, Si]] + 2[So, [Si, Ci]] - 2[Co, [Ci, Si]] 

- 2[Si, [So, Co]] + 2[Ci, [Co, So]] + 4[So, [Si, Co]] - 4[Co, [Ci, Sq]] 
+ [So, [So, Ci]] - [Co, [Co, Si]] - [So, [So, Co]] + [Co, [Co, So]] 

- [So, [So, Si]] + [Co, [Co, Ci]] - [Si, [So, Si]] + [Ci, [Co, Ci]]}. (10.8) 

Our final explicit example is E^-^\ obtained from 

+ S^'^eSI + S^'^Egl = ig[6u;^,^, (10.9) 
(Recall that S'^^^ = 0.) It has fewer terms than might have been anticipated: 

^Sn = '-^{[Bi, [Bo, [Bo, Si]]] - [Ci, [Co, [Co, Ci]]] + 2[So, [Co, [So, Si]]] 

- 2[Co, [So, [Co, Ci]]] - 2[Si, [Co, [Sq, Si]]] + 2[Co, [Si, [So, Si]]] 
+ 2[Ci, [So, [Co, Ci]]] - 2[So, [Ci, [Co, Ci]]] + 2[Si, [Sq, [Ci, Si]]] 

- 2[Ci, [Co, [Si,Ci]]] - 2[So, [Si, [So,Ci]]] + 2[Co, [Ci, [Co, Si]]] 

- [Si, [Ci, [Si, Ci]]] - [So, [Co, [So, Co]]] - 2[So, [Ci, [Si, Ci]]] 
+ 2[Ci, [So, [Si, Ci]]] + 2[Co, [Si, [Ci, Si]]] - 2[Si, [Co, [Ci, Si]]] 

- 2[So, [Co, [So, Ci]]] + 2[Co, [So, [Co, Si]]] + [Si, [Co, [Si, Co]]] 

- [Ci, [So, [Ci, So]]] - 2[So, [Ci, [Si, Co]]] + 2[Ci, [So, [Si, Co]]] 

- 2[Ci, [Co, [Si, So]]] + 2[Co, [Ci, [Si, Sq]]] - 2[So, [Co, [Si, Ci]]] 

+ 2[Co,[So,[Ci,Si]]]} (10.10) 

Although we can compute 

E(N) 

to any required order A^, we have not yet discovered 
a general iterative formula for E^^K A notable feature is the appearance in of the 

term 

(igh)''^ J2[A^, [A^, [..., [A, m] . . .]]], (10.11) 



where {Ai} = {So, Si, Co, Ci}, which, for example, results in a significant simplification 
of (10.8), and explains the absence of many terms in (10.10), because B3 = 0. The 
calculational methods so far developed for extracting information about matrix elements 
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require the knowledge of only the first few orders in any case, so the first few E^^^ should 
be sufficient for at least the initial stages of the finite- clement solution. 

The generalization to four dimensions is nontrivial but straightforward. Of course 
there are more Yang-Mills field components, and new structures emerge. The lowest-order 
interaction terms are given in Ref. 30. 

B. Construction of Yang-Mills Equations 

We finally must construct the lattice analogue of the continuum Yang-Mills equations, 
(9.9). Upon a moment's refiection, however, it is clear that the structure of the interactions 
in the Dirac equation carries over to the Yang-Mills equations, with multiplication by 
powers of the potentials replaced by nested commutators. The reason for this is as follows. 
We must solve the covariance equation 

SErn,n = ig[SuJm,n, Em,n] (10.12) 

for SF and 6E, where F and G are given by (9.16). Because necessarily the boson fields 
must be periodic, to do this we require that the number of spatial lattice sites, M, be odd. 
Then we find 

m+M 

SFm,n=ig[^U;m,n,Fm,n] + ^-^ J2 (-1)™^™' [<^^m',n - (^a;^'-l,n, i^m',n] , (10.13) 

m'=m+l 

and 

n 

SGm,n = ig[SuJm,n, Gm,n] - (-1)'^+'^ [5L0m,n' - 5Wm,n'-l , G^.n'] • (10.14) 

n'=0 

These have just the form as the transformation of 9 and 0, (9.19) and (9.39), respectively. 
Here, as in (9.37), we have adopted the initial conditions, at fixed spatial coordinate m, 

5Go = '-^[Sooo + 5u_i, Go]. (10.15) 

Without more ado, we transcribe the form of the interaction terms. For F we suppress 
the local variable n and write, for the term to be added to the left-hand side of (9.15a), 

^ m+M N / . , 

^^r^^ = l E (-1)-+-' Yl -^^ l^'"'- 1^'"'- ^ ■ • n'"^,^'"^] ■ ■ •]]]; (10-16) 

m'=m+l k=l , , ^ 

fc nested commutators 

and for G we suppress the local variable m and write, for the term to be added to the 
left-hand side of (9.15b), 

n N , . 



n' = l fe=l 

k nested commutators 

+ ii9hr^i-ir [Bo, [Bo, [. . .,[Bo, 4°)] . . .]]] . (10.17) 

N nested commutators 
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Here Im^ = —2Fm/h and ]Cn^ = —2Gn/h. In Ref. 29 we prove the required covariance 
statements: 



AT 



and 



fe=0 
N 



k=0 



The Yang-Mills equations are given by (9.15), (10.16), and (10.17): 

1 



h 



(-^m+l,n -^m,n) ~l" -^m,n — J 



m.ni 



h 



where 



iGm,n+l Gm,n) ~l~ ^m,n — 3m,ni 

00 

-^m,n — / ^ -^m,ni 
N=l 

00 



(10.18) 
(10.19) 

(10.20a) 
(10.20b) 

(10.21a) 
(10.21b) 



Ar=i 



Here, we easily derive from (10.16) and (10.17) first 

^ m+M 

= 9 E (-1)"^+"^' [e-*^'^^'"' (X^^} + Xm) e^^'^^-' - X^^, - X^\ , (10.22a) 



m'=m+l 



and 



n 

/C„ = - ^ (_l)n+n' Jg-i5/.B„, j^^(0) ^ ^^^^ _ ^(0) _ 



n'=l 



+ (-iy 



JghBo/2^{0)^-ighBo/2 _ j^{0) 



and then the difference equations 



and 



T e'^ghCm I „ighCm'j- , r„isfeC„i 771 ] 



h 



(10.22b) 



(10.23a) 



(10.23b) 



We are currently studying the solutions of these equation in two and four space-time 
dimensions. 
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